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Abstract 

Background: In both beef and dairy cattle, the majority of early embryo loss occurs within the first 14 days following 
insemination. During this time-period, embryos are completely dependent on their maternal uterine environment for 
development, growth and ultimately survival, therefore an optimum uterine environment is critical to their survival. The 
objective of this study was to investigate whether differences in endometrial gene expression during the mid-luteal 
phase of the estrous cycle exist between crossbred beef heifers ranked as either high (HF) or low fertility (LF) (following 
four rounds of artificial insemination (Al)) using the Affymetrix® 23 K Bovine Gene Chip. 

Results: Conception rates for each of the four rounds of Al were within a normal range: 70-73.3%. Microarray analysis 
of endometrial tissue collected on day 7 of the estrous cycle detected 419 differentially expressed genes (DEC) 
between HF (n = 6) and LF (n = 6) animals. The main gene pathways affected were, cellular growth and proliferation, 
angiogenesis, lipid metabolism, cellular and tissue morphology and development, inflammation and metabolic 
exchange. DEC included, FST, SLC45A2, MMP19, FADS! and GALNT6. 

Conclusions: This study highlights, some of the molecular mechanisms potentially controlling uterine endometrial 
function during the mid-luteal phase of the estrous cycle, which may contribute to uterine endometrial mediated 
impaired fertility in cattle. Differentially expressed genes are potential candidate genes for the identification of genetic 
variation influencing cow fertility, which may be incorporated into future breeding programmes. 



Background 

The failure of breeding females to become pregnant, in 
both dairy and beef cattle production systems, directly 
impacts the economic viability of these enterprises, and 
ultimately hinders genetic progress. Significant decreases 
in dairy cow fertility, ranging from 0.45% to 1% per annum, 
have been reported in cattle populations across the globe 
[1-3]. Following insemination the greatest increment of 
cow reproductive wastage occurs in the form of early em- 
bryo mortality with approximately 80% of this occurring 
within 14-16 days [4-6]. More specifically, previous studies 
have highlighted that the majority of early embryo loss 
typically commences around the mid-luteal phase of an 
estrous cycle i.e. day 7 of pregnancy [7,8] concurrent with 
the critical blastulation stage of embryo development [9]. 
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There is evidence of repeatable differences between cows 
in their ability to become pregnant. McMillan [10] reported 
a 65% difference in pregnancy rate at 60 days of gestation, 
following 6 consecutive in vitro embryo transfer events, be- 
tween two groups of cows. Differences in follicle wave dy- 
namics, duration of estrus, site of ovulation, or subsequent 
progesterone profiles were not found to contribute to the 
observed difference in pregnancy rate. Indeed, the authors 
suggested that "uterine" rather than "ovarian" factors may 
be responsible for the variation observed. This uterine ef- 
fect was also hypothesized in similar studies examining 
phenotypic differences between high and low fertility ani- 
mals [11-13]. Furthermore, data from our laboratory sug- 
gest a repeatability estimate of 0.18 for embryo survival in 
beef heifers [13] and heritability estimates for conception 
rate have been reported to exceed 0.20 [14,15]. 

The prerequisites to the establishment and mainten- 
ance of a successful pregnancy include a viable embryo, 
an appropriate steroidal environment and an optimally 
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functioning and receptive endometrium [16-18]. The endo- 
metrium plays a pivotal role in orchestrating the events 
that lead to fertilization, implantation and pregnancy. 
Throughout the estrous cycle and pregnancy, the endomet- 
rium is subjected to a host of functional and morphological 
changes, regulated by the hormones progesterone, estradiol 
and oxytocin [19]. The endometrium also functions to 
secrete a multitude of growth factors, proteins and cyto- 
kines, all of which constitute the histotroph, an import- 
ant source of energy and nutrition to a growing embryo 
in vivo [20-22]. 

Using conventional candidate approaches many stud- 
ies have examined bovine endometrial gene expression 
under various conditions; during early pregnancy in ani- 
mals that produced viable and non-viable embryos [23], 
in pregnant and cycling animals with artificially induced 
high, and normal systemic progesterone concentrations 
[24-26] and during the various phases of the estrous cycle 
[27]. Furthermore global endometrial gene expression ana- 
lyses have been conducted and include comparisons be- 
tween cycling and pregnant animals [28,29], fertile and sub- 
fertile animal strains [30,31], progesterone supplementation 
treatments [32], and specific estrous cycle phases [33,34]. 
Despite these efforts, endometrial gene expression of ani- 
mals characterized as either high or low fertility has not 
been investigated. Given the critical importance of day 7 
[7,8], we hypothesise that uterine endometrial gene expres- 
sion patterns will be different between high and low fertility 
heifers on day 7 of the estrous cycle. Thus, the objective of 
this study was to characterize differential gene expression 
profiles in endometrial tissue harvested on day 7 of the es- 
trous cycle from heifers ranked as either HF or LF fertility 
based on four successive inseminations and pregnancy 
diagnoses. Intercaruncular endometrial tissue was exam- 
ined due to the fact that caruncular endometrium lacks 
uterine glands which are essential to the exchange, trans- 
port and secretion of pertinent metabolites which consti- 
tute the uterine histoptroph and are required to support 
pregnancy [35,36]. 

Methods 

Ethics statement 

All experimental procedures involving heifers were li- 
censed by the Department of Health and Children, Ireland 
(licence number BlOO/846). Protocols were in accordance 
with the Cruelty to Animals Act (Ireland 1876, as amended 
by European Communities regulations 2002 and 2005) 
and the European Community Directive 86/609/EC and 
were sanctioned by the Institutional Animal Research 
Ethics Committee. 

Animal model 

Estrous cycles of reproductively normal nuUiparous cross- 
bred beef heifers {Bos taurus n = 120) were synchronized 



using two intramuscular administrations of 500 [ig 
of the prostaglandin analogue (PG), cloprostenol 
(Estrumate*, Schering-Plough Ltd., Shire Park, Welwyn 
Garden City, Hertfordshire, UK). Animals were visually 
observed for signs of estrous activity 3- to 5- times daily 
as described by Lynch et al. [7]. Only heifers observed to 
be in standing estrus were inseminated 6-18 hrs after on- 
set of heat [37]. Inseminations were carried out artificially 
by one trained technician. Heifers were given a single 
insemination of frozen-thawed semen, collected from a 
single ejaculate of one high fertility bull. Sire breed was 
Limousin and named Bolide (FL17). At the time of the 1st 
insemination, heifers were on average 20 months of age 
and weighed 440 ± 9.0 kg (Mean ± SEM). 

Using an Aloka SSD-500 V ultrasound scanner, fitted 
with a 7.5 MHZ transducer (Aloka Co. Ltd., Tokyo, Japan), 
pregnancy was diagnosed 28 days after insemination using 
the criteria set out by Kastelic et al. [38]. Following diag- 
nosis, all pregnant heifers received PG on day 28 to induce 
embryo loss. Six weeks after induced embryo loss all 
heifers were subjected to estrous reprogramming using a 
two-injection PG-regimen (11 days apart), inseminated 
and pregnancy scanned as described above. 

For the purpose of establishing an accurate high versus 
low heifer fertility model, this schedule was followed for 
a further two occasions. Thus, following four insemina- 
tions, animals that established a pregnancy on all four 
occasions were categorized as "HF" heifers while those 
achieving pregnancy on only one occasion were catego- 
rized as "LF" heifers. To eliminate the possibility of a phys- 
ical or anatomical abnormality that may have impeded 
heifers from becoming pregnant, animals with zero re- 
corded pregnancies were omitted from the study. 

After the fourth insemination, and subsequent pregnancy 
diagnosis, pregnant heifers were returned to estrous. Ap- 
proximately three months later, estrous cycles of animals 
were synchronized again in preparation for endometrial 
harvesting on D7. Figure 1 illustrates the timeline of events 
during the experimental period. 

Throughout the experimental period, animals had ad 
libitum access to grass silage supplemented with 2 kg of 
concentrates per heifer per day. Heifers were housed on 
concrete slats in groups of 15, at 2.5 m^ per heifer, for 
the duration of the study (15 months). Slaughter live- 
weight averaged 625 kg, with ECS of 4.0. Heifers were 
gaining weight during the course of inseminations in the 
region of 0.60 kg/day. 

Tissue sampling 

Animals from HF (n = 6) and LF groups (n = 6) were 
slaughtered on D7 in a licensed abattoir (KEPAK, 
Athleague, Co. Roscommon, Ireland). Following slaughter 
the reproductive tract and ovaries were checked for gross 
abnormalities but none were recorded. Uterine tissues 
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were opened longitudinally along the mesenteric border. 
Intercaruncular endometrial cross-sections approximately 
4 sq cm, and weighing 2.5 g, were harvested from the ani- 
mals by peeling from the underlying uterine myometrium 
from the middle-third of the uterine horn ipsUateral to the 
corpus luteum (CL) within 20 min of slaughter. 

Samples were washed in sterile PBS, and stored in 
RNAlatef at 4°C for 24 h before being transferred for long- 
term storage at -20°C. All surgical instruments used for tis- 
sue collection were sterilized and treated with RNA Zap 
(Ambion, Applera Ireland, Dublin, Ireland). In addition, on 
the day of slaughter CL diameter for each heifer was deter- 
mined using vernier calipers. 

Blood sampling 

Heifers were blood sampled via jugular venipuncture for 
subsequent measurement of progesterone at 0900 and 
2100 h commencing 24 h after PG for a cycle length. All 
blood samples were collected into 10 ml ethylenediamine 
tetraacetic acid (EDTA) heparinized Vacutainers (Becton 
Dickson Vacutainer Systems, Plymouth, UK). Samples were 
held in iced water untO centrifuged at 1500 x ^ at 4°C for 
15 mins after which plasma was extracted and stored in 
sterile 7 ml vials at -20°C until assayed. 

Progesterone assays 

Progesterone profiles for each of the six heifers within 
HF and LF groups were established. Concentration of 
progesterone was measured in plasma as the mean of 
the two samples taken on each cycle day of the previous 
cycle and on 7 days prior to slaughter using the Coat-a- 
Count assay procedure (Coat-a-Count Diagnostic Products 
Corporation, Los Angeles, CA, USA) with each sample 
tested in duplicate. The inter-assay and intra-assay coeffi- 
cients of variation for low, medium and high control 



samples were 17.4% and 4.4%, 5.6% and 28.4%, and 4.2% 
and 4.9% with mean concentrations of 0.24, 2.54 and 
7.21 ng/mL, respectively. The minimum detectable limit 
for this assay was 0.06 ng/mL. 

RNA extraction and quality analysis 

Total RNA was prepared from 100-200 mg of endometrial 
tissue using the TRIzol reagent (Sigma-Aldrich Ireland 
Ltd., Dublin, Ireland). Tissue samples were homogenized 
in 3 ml of TRIzol reagent and chloroform, and subse- 
quently precipitated using isopropanol (Sigma-Aldrich 
Ireland Ltd., Dublin, Ireland). RNA samples were stored 
at -80°C. Samples of RNA, (20 |ig), were purified and 
treated for contaminating genomic DNA using RNeasy 
clean-up kits in accordance with manufacturer's guide- 
lines supplied (QIAGEN, Crawley, West Sussex, UK). 
This protocol included an on-column DNase treatment 
step. RNA quality and quantity were assessed using auto- 
mated capillary gel electrophoresis on a Bioanalyzer 2100 
with RNA 6000 Nano Lab-chips according to manufac- 
turer's instructions (Agilent Technologies Ireland, Dublin, 
Ireland). Absorbance ratios (28S/18S) and RNA integrity 
values recorded for all RNA samples extracted post 
clean-up ranged between 1.8 and 2.0, and 7.5 and 9.8, 
respectively. 

MIcroarray hybridization 

Gene expression was determined using a 24,027 probe 
set bovine oligonucleotide array (Affymetrix"), represent- 
ing -23,000 bovine transcripts based on the original 
mapping using Unigene build 57 (March 24, 2004). RNA 
from each heifer was hybridized to a separate array. All 12 
RNA samples were hybridized and scanned by the German 
Resource Centre for Genomic Research (RZPD), Germany, 
according to the manufacturer's instructions. 
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Microarray analysis 

All microarray analyses including preprocessing, normal- 
ization and statistical analysis were carried out using R 
(R, 2007) version 2.6 and Bioconductor [39] version 2.1 as 
previously described by [40]. Data were quality assessed 
before and after normalization using a number of in-built 
quality control methods implemented in the Bioconductor 
affycoretools and associated packages to identify problems 
if they existed with array hybridization, RNA degradation 
and data normalization. Microarray data were prepro- 
cessed using the mmgMOS normalization method [41,42] 
using the default settings and differential expression (DE) 
was calculated using the pumaDE method both imple- 
mented in the Bioconductor package "puma" [42-45] . The 
puma method uses a Bayesian hierarchical model to 
calculate the probability of positive likelihood ratio 
(PPLR). The PPLR associates probability values of genes 
being differentially expressed, which is a measure of false 
positive detection of DE, to each ratio and generates lists 
of genes ranked by the probability of DE. This PPLR 
statistic was converted into "P-like values" using the 
recommended formula in the puma method prior to sub- 
sequent analysis. 

As many of the original annotations for the Affymetrix 
bovine chip are erroneous [6,46], remapped annotations 
were determined using the "bovinedaiplusv6cdf' chip def- 
inition file (CDF). This annotation is based on the CDF- 
Merger procedure as described by De Leeuw et al. [47], 
which generates a hybrid CDF based on the standard Affy- 
metrix CDF (version 26) and the custom Brainarray (ver- 
sion 11.0.1) CDF. This re-mapped annotation includes 
mapping to all RefSeq (mature RNA protein coding 
transcripts and vaUdated complete coding sequences in 
GenBank). Annotations were also supplemented by inter- 
rogating the Ensembl Bos taurus database version 46 using 
the BioMart package in Bioconductor and manual annota- 
tion where possible with recent entries in Entrez Gene. 

Pathway analysis 

To examine the molecular functions and genetic net- 
works, the microarray data were further analyzed using 
Ingenuity Pathway Analysis (v. 8.8, Ingenuity Systems, 
Mountain View, CA; http://www.ingenuity.com), a web- 
based software application that enables identification of 
over-represented biological mechanisms, pathways, and 
functions most relevant to experimental datasets or genes 
of interest [40,48-50]. 

A dataset containing gene identifiers and corresponding 
expression and P-like values was uploaded into IPA. Briefly, 
each identifier was mapped to its corresponding gene 
object in the Ingenuity knowledge base. A P-like value 
oi P < 0.05 from the puma analysis was set to identify 
genes whose expression was statistically significantly 
up- or down-regulated. These genes, called "focus" genes. 



were overlaid onto a global molecular network developed 
from information contained within the Ingenuity know- 
ledge base. Networks of these focus genes were then al- 
gorithmically generated based on their connectivity. 
Network analysis returns a score that ranks networks 
according to their degree of relevance to the network 
eligible molecules in the dataset. The score takes into 
account the number of network eligible molecules in 
the network and its size, as well as the total number of 
network eligible molecules analyzed and the total num- 
ber of molecules in the knowledge base that could po- 
tentially be included in networks. 

RT-qPCR analysis 

The microarray results were validated by carrying out 
RT-qPCR on 18 genes. Candidate genes were chosen 
based on the following criteria; those that were top rank- 
ing in our microarray DEG list, genes with known func- 
tional importance in uterine mediated sub-fertility which 
were either up- or down-regulated and genes which 
were not differentially expressed between the two treat- 
ment groups. 

Using the same RNA samples that were analyzed in 
the microarray studies, first strand cDNA was synthe- 
sized using the High Capacity cDNA Reverse Transcrip- 
tion kit according to manufacturer's instructions (Applied 
Biosciences, Ireland). Purified total RNA (1 |.ig) was re- 
verse transcribed using random hexamers. The converted 
cDNA was quantified by absorbance at 260 nm, diluted to 
50 ng/fil working stocks and stored at -20°C, for subse- 
quent analyses. 

Analysis of putative reference genes for RT-qPCR stud- 
ies was carried out using GeNorm version 3.5 Microsoft 
Excel Add in (Microsoft, Redmond, WA) [51]. The stabil- 
ity of the expression of several cited reference genes in- 
cluding, ribosomal protein L15 [52], 18 s ribosomal RNA 
[53], ubiquitin [54], glyceraldehyde phosphate dehydro- 
genase and |3-actin [55,56], was investigated across all 
samples in this study. Similar to Coyne et al. [54], ubiqui- 
tin (at an optimal concentration of 2.5 i^M) exhibited the 
greatest stability during qPCR analysis of endometrial 
mRNA samples analyzed, with an M value of 0.022. Based 
on a recommended cut-off lvalue of 0.15; ubiquitin was 
selected as a single standard reference gene for these ex- 
periments as the use of additional reference genes did not 
contribute to a more accurate normalization factor. 

Primers were designed, to span exon boundaries where 
possible, using the PrimerS software programme [57] 
and oligos were aligned by Basic Local Alignment Search 
Tool (BLASTN) on the National Centre for Biotechnol- 
ogy Information (NCBI) web page, to verify their identity 
and homology to the bovine genome (http://www.ncbi. 
nlm.nih.gov/BLAST/). All oligonucleotides were commer- 
cially synthesized as highly purified salt-free products by 
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Sigma Aldrich Ireland Ltd. Primers were first tested using 
end point PCR to optimize amplification conditions. All 
amplified PCR products generated in this study were 
purified using the PCR purification kit (Roche, Basel, 
Switzerland) and sequenced (Macrogen; Nucleics Pty 
Ltd, Bendigo, Australia) to verify their identity. Primer 
sequences used in this study are listed in Table 1. 

Primer concentrations were optimized for each gene 
by titrating 5, 10, and 20 |iM per primer. The most suit- 
able primer concentration was chosen based on four cri- 
teria in order of decreasing importance: i) a clear distinct 
melt curve absent of any additional peak(s) caused by 
non-specific binding, ii) a curve within the temperature 
range 75-85°C, iii) the primer concentration producing 
the lowest threshold cycle number (Ct) and lastly, iv) repli- 
cation amongst Ct values and melting temperatures (Tm). 
Subsequently, efficiencies of chosen primer concentrations 
were determined over a 5-fold dilution series, whereby 
cDNA was diluted into working solutions: stock, 1:2, 1:4, 
1:8, 1:16, and RT-qPCR assays carried out. This was re- 
peated for every gene. The and amplification efficiency 
(£) values for RT-qPCR were calculated from linear re- 
gression analysis of log (input cDNA) versus Ct plot. The 
slope for each set of standards was used to determine E = 
^g(-i/siope) _ ^ Slopes, amplification efficiencies and R^ 
estimates for individual genes are reported in Table 2. 
Only primers with PCR efficiencies between 90% and 
110% were used. 

Each RT-qPCR reaction was carried out in a 96-well 
plate format with a total volume of 20 |il, containing 
1 \x\ cDNA, (10 ng/\A), 10 ^1 Fast SYBR» Green Master 
Mix (Applied Biosystems, Ireland), 1 i^l forward and reverse 
primers and 8 |il nuclease-free H2O. Performance of RT- 
qPCR was carried out using the Applied Biosystems Fast 
7500 v2.0.1 with the following cycling parameters: 95°C for 
10 min followed by 40 cycles of 95°C for 15 s and 60°C for 
60 s, followed by amplicon dissociation (95°C for 15 s, 60°C 
for 60 s, 95°C for 15 s and 60°C for 15 s). Dissociation 
curves were examined for the presence of a single PCR 
product. The software package GenEx 5.2.1.3 (MultiD Ana- 
lyses AB, Gothenburg, Sweden) was used for efficiency 
correction of the raw cycle threshold (Ct) values, inter- 
plate calibration based on a calibrator sample included 
on all plates, averaging of replicates, normalization to 
the reference gene and the calculation of quantities 
relative to the greatest Ct. Expression of each target 
gene was normalised to the reference gene and relative 
differences in gene expression were calculated using the 
2'^'^c'r method [58]. 

Statistical analysis 

All data were analyzed using the Statistical Analysis Sys- 
tems software package (SAS Inst. Inc., Cary, NC) version 
9.1. Data from RT-qPCR studies were tested for adherence 



to normality using PROC UNIVARIATE (SAS, 2003). 
Non-normal data were subsequendy transformed using 
the best fit fimction as described by PROC TRANSREG 
(SAS, 2003). Differences in mean values between the two 
groups (HF and LF) were tested using ANOVA (PROC 
MIXED). Animal within treatment was used as the error 
term. The Tukey critical difference test was used to deter- 
mine statistical difference between LF and HF mean values. 
The CORR procedure of SAS (PROC CORR, SAS 2003) was 
used to determine correlations between microarray and RT- 
qPCR data. Pearson correlation coefficients were estimated 
for each individual gene across all animals (n = 12). A P value 
of P < 0.05 was considered to be statistically significant. Data 
collected from CL diameter measurements were tested for 
adherence to normality using PROC UNIVARLVTE (SAS, 
2003). CL differences in mean values between the two 
groups (HF and LF) were tested using ANOVA (PROC 
MIXED). Animal within treatment was used as the error 
term. For the analysis of progesterone profiles individual 
profiles were normalized relative to day of estrus (Day 0). 
The effect of fertility status "HF" versus "LF" was established 
using a repeated measured analysis (PROC MIXED; SAS). 

Results 

Animal model 

Embryo survival rates were 73.3%, 71.7%, 73.3% and 
70.0% for A.I. rounds 1, 2, 3 and 4 respectively. A total 
of 31 heifers qualified as HF or LF; 15 HF and 16 LF, of 
which three of these were eliminated from the study due 
to the presence of ovarian abnormalities detected at 
ultrasound scanning. Pregnancy rate for LF heifers was 
consistent across all four replicates. Six HF and 6 LF heifers 
were randomly chosen within their respective fertility 
groups for slaughter on D7. The mean inter-estrous inter- 
vals in a previous recorded estrous cycle were 20.17 ± 0.96 
and 20.83 ± 0.96 days (P > 0.10) for the HF and LF heifers, 
respectively. At day of slaughter, mean CL diameters were 
22.58 ± 3.48 (SD) mm and 23.55 ± 4.4 (SD) mm for HF and 
LF heifers, respectively, i.e., there was no significant differ- 
ence in CL diameter between fertility groups {P > 0.10). 

Progesterone profiles 

There was no effect of fertility status, or interaction efiect 
of fertility status and day of cycle (P > 0.10), on the concen- 
tration of progesterone. On the day of slaughter plasma 
concentrations did not differ between the high and low fer- 
tility groups (HF 5.96 ng ml ^ LF 5.65 ng mV^, P = 0.589). 

Microarray differential gene expression 

A total of 419 genes were found to be differentially 
expressed between LF and HF (n = 6 vs. 6). Of these, 171 
were up-regulated and 248 down-regulated in the LF com- 
pared with HF heifers, respectively. Transcript abundance 
differences between LF and HF groups resulted in fold 
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Table 1 Bovine specific oligonucleotide forward and reverse primer sequences (5-3') and PCR product length 

Gene name Sequence Accession number Amplicon size (bp) 



18S/285 


F: 5'- TGCTCTCGCAAACCrAACa-3' 
R: 5'- CACTAAGCACTCGCATCCA-3' 


DQ222453 


159 


ACTA2 


F: 5'- ACTGGGACGACATGGAAAAG -3' 
R: 5'- TACATGGCTGGGACATTGAA-3' 


BT021508 


166 


AaB 


F: 5'- ACTrGCGCAGAAAACGAGAT-3' 
R: 5'-CACOTCACCGTTCCAGTrr-3' 


BT030480 


121 


APEH 


F: 5'- CAAGAGCATGCGCAGTATGT -3' 
R: 5'- GTAGAGCTGCAAAGCCCATC-3' 


BCl 23400 


181 


CELAl 


F: 5'- GGAACCATCCTGGCTAACAA -3' 
R: 5'- CATGGTGGTCTTCACAGTGG -3' 


BC1 49525 


165 


DONSON 


F: 5'- TGTGTTGGTGAAGGGAATGA -3' 
R: 5'- AGAGGGTTGGTGGAAGTCCT -3' 


BCl 33573.1 


107 


EST 


F: 5'- TAAATGAGAGACCCGCCAAC-3' 
R: 5'- CCCCAG'nTCTGTCCTGTGT-3' 


AY775795 


171 


GALNT6 


F: 5'- GACCACGTCTTGGACCTCAT-3' 
R: 5'- AGCTCAGCTGGGGTGTAGAA-3' 


NM_001015534 


146 


GAPDH 


F: 5'- GGGTCATCATCTCTGCACCT-3' 
R: 5'- GGTCATAAGTCCCTCCACGA-3' 


NM_001 034034 


176 


GJAl 


F: 5'- - CAACATGGGTGACTGGAGTG 3' 
R: 5'- GCAGGATTCGGAAAATGAAA -3' 


BT021508 


110 


IL33 


F: 5'- 1 IGI 1 1 IGGAGGATGGAAGC -3' 
R: 5'- TTTGTGGGGCTCAGGTITAC -3' 


BCl 23562 


163 


MMP19 


F: 5'- TGGACGTrATCCCCTCAGTC-3' 
R: 5'- GTCCATGGTTCATGCTTGTG -3' 


BCl 23722 


119 


M0SC2 


F: 5'- GCAGTGCTTTTGAGGAGGAC-3' 
R: 5'- GGATCACACAGGCGGTAACT-3' 


NM_001 076380 


169 


NMB 


F: 5'- ACATGACGACATGGCTGAAA-3' 
R: 5'- ACTrCAACAGGGAAGCGAGA-3' 


NM_001 075270 


185 


NPPC 


F: 5'- GAGGCAACAAGAAGGGTTTG -3' 
R: 5'- CTGATGACCAAGGGTGACCr -3' 


BCl 23399 


149 


PPARA 


F: 5'- rrGTGGCTGCrATCATTrGC-3' 
R: 5'- AGAGGAAGACGTCGTCAGGA-3' 


AF229356 


135 


RAB3B 


F: 5'- TGGGCGGAGATTCATTTTAC -3' 

n. D \JJr\rAr\A\JJ 1 u 1 0^r\ 1 OUVJJ 1 O 1 O J 


BCl 12795 


144 


RPU5 


F: 5'- TGCATAAGCACAGGGAGATG-3' 
R: 5'-CTGGAGAGTATTGCGCCTTC-3' 


BT020706 


134 


SERPl 


F: 5'- GTCCCTCTGGGTGAATCTGA-3' 
R: 5'- TCACTAATTGCCAGGGGTrC-3' 


NM_1 74460 


158 


SLC1A3 


F: 5'-CATCCATGCTGTCATTGTCC-3' 
R: 5'-ATCTGGTAACGCGTTTGTCC -3' 


BCl 201 25 


188 


SLC45A2 


F: 5'- CATGCCCTCTTCACAGGTTT-3' 
R: 5'- AGTGGGGCTTCAGGGATACT -3' 


XM_001251343 


179 
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Table 1 Bovine specific oligonucleotide forward and reverse primer sequences (5-3') and PCR product length 

(Continued) 

TGFBlll F: 5'- CCTGCAATAAACCCATTGCT-3' NM_001035313 162 

R: 5'- AGAAGCGCTCGAAGTAGCAC-3' 
UBQ F: 5'- TAGAACAGTTGGTGGCGAAA-3' BC1 02888 121 

R: 5'-GAAGAGTGGGCrGAaGAGG-3' 

18S/28S, 18S/28S ribosomal RNA; ACTA2, actin alpha 2; ACTB, beta-actin; APEH, N-acylaminoacyl-peptide Hydrolase; CELAl, chymotrypsin-like elastase family mem- 
ber 1; D0N50N, downstream neighbor of SON; FST, follistatin; GALNT6, UDP-N-acetyl-alpha-D-galactosamine:polypeptide N-acetylgalactosaminyltransferase 6; 
GAPDH, glyceraldehyde-3-phosphate dehydrogenase; GJA1, gap junction protein, alpha 1; IL33, interleukin 33; MMP19, matrix metallopeptidase 19; fV10SC2, MOCO 
sulphurase C-terminal domain containing 2; NMB, neuromedin B; NPPC, natriuretic peptide C; PPARA, peroxisome proliferator-activated receptor alpha; RAB3B, 
RAB3B member ras oncogene family; RPL15, ribosomal protein LIS; SFRP1, secreted frizzled-related protein 1; SLC1A3, solute carrier family 1 member 3; SLC4SA2, 
solute carrier family 45 member 2; TGFBlll, transforming growth factor beta 1 induced transcript 1; UBQ, ubiquitin. 



changes ranging from 6.6-fold down to 8-fold up-regulated 
in LF animals. The microarray data have been deposited in 
NCBI's Gene Expression Omnibus [59] and are accessible 
through GEO Series accession number GSE29853. Hier- 
archical clustering of differentaiUy expressed genes is pre- 
sented as a heatmap and dendogram in Additional file 1: 
Figure SI. 

Pathway analysis 

Of the 419 DEG, a total of 227 genes were successfully 
mapped to a molecular/biological pathway and/or cat- 
egory in the IPA database, while 202 of these were net- 
work eligible using IPA. Among the mapped DEG, 73 
were up-regulated (Additional file 2: Table SI) and 154 
down regulated (Additional file 2: Table S2). 



Table 2 Efficiency variables for individual RT-qPCR genes 



Gene 


Optimum [primer] pIVI 


Slope (-) 




Efficiency 


AaA2 


5 


3.62 


0.99 


1.890 


APEH 


20 


3.18 


0.99 


2.064 


CEiAl 


10 


3.54 


0.99 


1.917 


D0N50N 


10 


320 


0.99 


2.054 


FST 


2.5 


3.60 


0.97 


1.896 


GALNT6 


5 


3.29 


0.99 


2.103 


GJAl 


20 


3.31 


0.98 


2.003 


IL33 


5 


3.33 


0.96 


1.997 


MMP19 


10 


3.52 


0.99 


1.925 


M0SC2 


10 


3.14 


0.99 


2.082 


NMB 


10 


3.21 


0.99 


2.049 


NPPC 


5 


3.63 


0.99 


1.885 


PPARA 


5 


3.33 


0.99 


1.997 


RAB3B 


20 


3.51 


0.99 


1.927 


SFRPl 


5 


3.70 


0.97 


1.863 


SLC1A3 


10 


3.46 


0.97 


1.946 


SLC45A2 


10 


3.54 


0.95 


1.916 


TGFBlll 


20 


3.29 


0.96 


2.013 


UBQ 


2.5 


3.20 


0.99 


2.054 



Biological functions 

Biological categories with the largest number of up regu- 
lated genes included DNA replication, recombination and 
repair, nucleic acid metabolism and carbohydrate metabol- 
ism. Categories with the largest number of down-regulated 
genes were organ morphology, and connective tissue de- 
velopment and function. Of the top 20 most statistically 
significantly over-represented biological categories, DNA 
replication, recombination and repair had the greatest ratio 
of up- to down-regulated genes (Figure 2). Pathways with 
the greatest number of DEG, including their respective 
number of DEG, were cellular growth and proliferation 
(n = 57), inflammatory disease (n = 55), cell death (n = 49), 
cellular development (n = 43), small molecule biochemis- 
try (n = 37), cellular morphology (n = 36) and tissue devel- 
opment (n = 36) as shown in Table 3. 

Canonical pathways 

Canonical signaling pathway analysis uncovered genes 
with functions in ILK-signaling, TR/RXR activation, regu- 
lation of actin based motility by Rho and Integrin signaling 
(Table 4). Genes associated with canonical signaling path- 
ways were down-regulated in LF animals for all statistically 
significant pathways mapped with the exception of TR/ 
RXR activation where the ratio of up- to down-regulated 
genes was uniform. Canonical metabolic pathways over- 
represented within the microarray data included fatty acid 
biosynthesis, o-glycan biosynthesis and purine metabolism. 
There were more genes up-regulated in canonical meta- 
bolic than canonical signaling pathways with the greatest 
ratio of up- to down-regulated genes expressed in the 
metabolic pathway: o-glycan biosynthesis (Table 4). 

Networks 

Using IPA a total of 19 gene networks were identified, 12 
of which had 13 to 25 focus genes among DEG (Additional 
file 2: Tables SI and S2). The 12 top networks are listed in 
Table 5. Lipid metabolism featured in three of the top 12 
networks. In addition, organ/tissue/cell morphology and 
development appeared a central biological theme over- 
represented among DEG. Illustrations of gene interac- 
tions among DEG contained within the top two scoring 
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Tissue Development 
Skelelal & Muscular System Development & Function 
Cell Morphology 
Tumor Morphology 
Small Molecule Biochemistry 
Nucleic Acid Metabolism 
Molecular Transport 
DNA Replication, Recombination, & Repair 
CeUular Development 
Hepatic System Development & Function 
Comiective Tissue Development & Function 
Organ Morphology 
Organ Development 
Skin Development and Funcdon 
Gene Expression 
Ceh DeaUi 

Dermatological Diseases & Conditions 
Cellular Growth & Proliferadon 
Carbohydrate Metabohsm 
Endocrine System Development & Function 





0 0.5 1 1.5 2 2.5 3 3.5 4 

-log (P-value) 

Figure 2 Classification of DEG according to top 20 molecular and cellular functions, most significantly affected by endometrial related 
sub-fertility, using IPA. The red/green bars indicate tlie likeliliood [-log (P-value)] tliat the specific molecular and cellular function category was 
affected by endometrial related sub-fertility compared with others represented in the list of DEG. The proportion of up- and down-regulated genes in 
each group is represented by the red and green segments on each bar, respectively. 



Table 3 Biological categories from IPA analysis with the largest number of DEG 



Biological category 



Gene symbol 



Cellular growth and proliferation 
nflammatory disease 
Cell death 

Cellular development 
Small molecule biochemistry 
Cell morphology 
Tissue development 



AaB, ADCY3, CAT. CD58, CIAOI. CN0T8, C0UA2, aSL2, CXCL2, aSL2, DAP, DCK, EMP3, ENPPl, ETFDH, FASN, 
FHL2, FSl FUS, GADD45B, GJAl, GLI3, HMGBl, IL65T, ILF3, IMPDH2 LGAL51, LGALS3, LTBPl, LY6E, MAGEDl, 
MAP3K7IP2, MEST, MMP19, NMB, ODCl, PARK?, PPARA, PHLDAl, PPARA, PRKRA, PRMTl, PSMBIO, PSMD2, PTPRK, 
RBBP7, SERPINHl SFRPl, SLC1A3, S0X6, TCFU, TCFBUl, T0B2, TP53in, UCP3, WFDCl, ZNF259 

ACJA2, ACFB, ADCY2, ADCY3, APEH, C1QTNF6, CACNB2, CAT, CD58, C0R02A, CXC12, DAP DCK, DSC2, EEF2, 
ENPPl, FAU, FNBPl, GADD45B, GALNT2, GU3, GLUL, GSDMB, HMGBl, IFI6, IL33, IL6ST IMPDH2, LGALSl, LGALS3, 
LPHN2, LTBPl, MYLK NSF ODCl, 0RC5L, PARP4, PCMTDl, PPARA, PRMTl, PTGIS, PTPRK RAB12, RARRES2, 
RGNEF RPS3, SCG5, SLC1A3, SLC25A24, SLC45A2, S0X6, SRPKZ TCFU, TNIK WFDCl 

ACTB, AaCl, BAG3 CACNB2, CAT CJSL2, CXCL2, DAR EMP3 FASN, FAU, FHL2, EST FUS, GADD45B, GALNT2, GIMAP5, 
GJAl, GLI3, HMGBl, HSPBl, IF16, IL6ST LGALSl, LGALS3, LTBPl, MAGEDl, MAP3K71P2, MYLK NPPC NSF ODCl, PARK7, 
PARP4, PHLDAl, PPARA, PRKRA, PTGIS, OKI, RCAN2, RPS3, SCG5, SFRPl, SLC1A3 SLC25A24 S0X6, SRPK2, TCFU, TGFBlll 

ARHGDIG, ARHGEFl 1, BOQ CAT CJSL2, CXCL2, ENPPl, FHL2, FHL3, FLNC EST FUS GADD45B, GEFT, G1MAP5, 
GJAl, GL13, HMGBl, IL33, IL6ST LGALSl, LGALS3, LTBPl, MAGEDl, MARCKS, MYLK NPPC NROBl, ODCl, PPARA, 
OKI, RARRES2, RCAN2, RNF128, SFRPl, SLC1A3, SM0C2 S0X6, STIMl, TCFU, TGFBlll, TNIK T0B2 

ACATl, ANKRD26, CAT CMPKl, DCK ERH, ETFDH, FASN, EST GALNT2, GJAl, GLUL, HMGBl, IL6ST LGALSl, LRAT 
MARCKS, NMB NPPC ODCl, PAICS, PARK7, PCCB, PPARA, PRMTl, PTGIS OKI, RAB3B SCG5, SLC1A3, SLC25AU, 
S0X6, SRD5AI, SULTIAI, SV2A, TGFBlll, UCP3 

ACTA2, ACTCl, ANXA6, ARHGDIG ARHGEFl I, CTSL2, CXCL2, DPYSL3, ENPPl, FASN, FERMT2, FHL3, FLNC EST 
GJAl, HMGBl, HSPBl, IL6ST LGALSl, LGALS3, LRAT MARCKS, NPPC NTO4, ODCl, PHLDAl, PPARA, RGNEF, 
RNFI28 SCG5, SERPINHl, SFRPl, SLCIA3, S0X6, TGFBlll, TNIK 

ACTA2, ACFCl, ALXl, ARHGDIG, BOC CAT CD58 C0L16A1, CXCL2, DSC2, ENPPl, FASN FHL3, GADD45B, GEFT, 
GJAl, GL13, HMGBl, IL33, IL6ST LGALS3, LY6E MARCKS, MEST MMP19, NPPC NROBl, NTN4, PHLDAl, PPARA, 
PTPRK SEMA5A, SFRPl, S0X6, TCFU, TGFBlll 
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Table 4 Enriched canonical pathways in endometrial mRNA from HF and LF heifers 

Pathways Genes 



% DEG 



P-value 



Canonical signalling 

Regulation of Actin-based Motility by Rho 
LK Signalling 

Caveolar-mediated Endocytosis Signalling 
Polyamine Regulation in Colon Cancer 
Cellular Effects of Sildenafil 
RhoA Signalling 

Mechanisms of Viral Exit from Host Cells 

TR/RXR Activation 

Virus Entry via Endocytic Pathways 

Semaphorin Signalling in Neurons 

HMGB1 Signalling 

Thrombin Signalling 

htegrin Signalling 

Sphingosine-1 -phosphate Signalling 
Germ Cell-Sertoli Cell Junction Signalling 
CXCR4 Signalling 
Canonical metabolic 

Purine Metabolism 

Glycan Biosynthesis 

Fatty Acid Biosynthesis 

Sulphur Metabolism 

Pantothenate and CoA Biosynthesis 

Pyrimidine Metabolism 



MYLK, ACTB, ACTA2, RHOU, ACTCl, FNBP 

TGFBin, FLNQ ACrN2, ACTB, AaA2, FERMT2, RHOU, AaCl, FNBPl 

FLNC ACTB ACTA2 ACTCl, C0PB2 

PSMD11, PSMBIO, PSMD2, ODCl 

MYLK, ADCY2, ACTB ACTAZ ADCY3, ACTCl 

MYLK, ACTB ACTA2, ARHGEFl 1, ACTCl 

ACTB ACTAZ ACTCl 

UCP3, RAB3B, FASN, RCAN2 

FLNC ACTB ACTAZ ACTCl 

DPYSL3, RHOU, FNBPl 

HMGBl (includes EQ3146), RBBP7, RHOU, FNBPl 
MYLK ADCYZ ADCY3, RHOU ARHGEFl 1, FNBPl 
MYLK ACFB AaAZ RHOU, ACTCl, FNBPl 
ADCYZ ADCY3, RHOU FNBPl 
ACTB ACTAZ RHOU ACTCl, FNBPl 
ADCYZ ADCY3, RHOU ARHGEFl 1, FNBPl 



NSF, ADCYZ ENPPl, IMPDHZ DCK, ATP 13 AS, ADCY3, 
PDE6C PAICS, P0LR2H, ACTCl 



GALNT2, GYLTLIB, GALNT6 
TASK PCCB 
SULTlAl, 5U0X 
ENPPl, DPYSL3 

ENPPl, DCK, DPYSL3, POLRZH, CMPKl 



6.5 
4.3 
6.0 
7.0 
3.9 
4.5 
6.8 
4.1 
4.2 
5.8 
4.1 
2.9 
3.0 
3.6 
3.2 
3.0 

2.5 

6.3 
3.9 
3.3 
3.1 

2.2 



0.0007 
0.0017 
0.0024 
0.0041 
0.0058 
0.0126 
0.0135 
0.0245 
0.0245 
0.0288 
0.0302 
0.0331 
0.0331 
0.0407 
0.0417 
0.0447 

0.0019 

0.0077 
0.0098 
0.0316 
0.0347 
0.0457 



Genes marked in bold are up-regulated. 



networks can be seen in Figures 3 and 4. Biological path- 
ways; lipid metabolism, cell growth and proliferation, and 
tissue development and function, were repeatedly featured 
pathways that constituted these top scoring networks. 

RT-qPCR analysis 

Eighteen genes were validated by real-time RT-qPCR 
(Table 1). There was moderate to good consistency be- 
tween methodologies for direction and magnitude of differ- 
ential gene expression among genes analyzed. Correlation 
coefficients exceeded 0.60 in fourteen of the eighteen genes 
validated (Figure 5, Additional fQe 2: Table S3). 

Discussion 

The animal model generated in this study, is the first of 
its kind. Two groups of heifers consistently divergent in 
conception rate; HF and LF were successfully generated 
and endometrial gene expression examined. We identi- 
fied key genes and pathways potentially contributing to 
endometrial related conception rate variance, the most 



extreme of which had no previously known involvement 
in endometrial function, including cellular growth and 
proliferation NPPC and GJAl; angiogenesis MMP19 and 
HMGBl; lipid metabolism FASN and PPARA; cellular and 
tissue morphology and development FST and TGFBIII; 
inflammation IL-33; and metabolic exchange SLC1A3 and 
SLC25A24. 

Several studies have highlighted the vital role progester- 
one plays in early embryo development to the extent that 
decreased conception rates were observed in heifers with a 
delayed postovulatory progesterone peak [60]. Furthermore, 
it has been well documented that progesterone influences 
endometrial and oviductal function [61,62]. In the present 
study, progesterone concentrations were within the normal 
range for both HF and LF heifers, and did not vary between 
groups. In addition, CL diameter measurements were not 
different between HF and LF animals and were consistent 
with observations from other studies examining CL diame- 
ters during this period of the estrous cycle [63] . The high 
conception rates achieved across successive breedings 
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Table 5 Networks generated from endometrial gene expression data of HF versus LF heifers by IPA 

Network ID Top functions Molecules In network 



Score 



Focus 
molecules 



Lipid Metabolism, Small Molecule Biochemistry 



Cellular Growth and Proliferation, Connective 
Tissue Development and Function, Skeletal and 
IVluscular System Development and Function 



Carbohydrate Metabolism, Haematological 
Disease, Metabolic Disease 



Cell Death, Gene Expression, Lipid Metabolism 



Cell Morphology, Inflammatory Response, 
Lipid Metabolism 



Drug Metabolism, Small Molecule Biochemistry, 
Cell-To-Cell Signalling and Interaction 



Embryonic Development, Organ Development, 
Organ Morphology 



Developmental Disorder, Neurological Disease, 
Cell Death 



Cellular Development, Cell Cycle, Connective 
Tissue Development and Function 



Post-Translational Modification, Cell Cycle, 
Gene Expression 



Cancer, Cardiovascular Disease, Cell-To-Cel 
Signalling and Interaction 



AaB, ADCY2, ADCY3, ANXA6, ARHGDIQ CACNB2, Calmodulin, 45 
CN0T8, CTSL2, ERK, F Actin, FAU, FKBP7, FKBPIO, FLNC FSH, 
GJAI, GUCY, hCG, IFI6, IBS, Lh, MARCKS (includes EG:4082), 
NPPC NROBI, NTN4, Peptidylprolyl isomerase, PHLDAI, Pkc(s), 
PPIH, Rock, RPLPI, RPLP2, SRD5A1, TP53II / 

AaA2 Alpha Actinin, CAT, CIAOI, C0UA2, Collagen type I, 45 
Collagen(s), C0PB2, CSRPI, CSRP2, CXCL2, ENPPI, FERMT2, 
FHL2 FHL3, EST, G-Actin, GADD45B GEFT, HSPBI, Integrin, 
Laminin, LTBPI, MAGEDI, MMP19, MYLK, NFkB (complex), Pak, 
Pdgf, PJAI, RARRES2, RITI, SERPINH1, Tgf beta, TGFB1II 

26 s Proteasome, Akt, AMPK, EIF2C4, EL0VL5, FA5N, GLUL, 44 
HISTONE Histone hS, Histone h4, HMG20B, Hsp90, MAP3K7IP2, 
MEDI3, MED27, N-cor, ODCI, PPARA, PSMBIO, PSMD2, 
PSMDII, PTPRK, RAB3B RBBP7, RCAN2, RPL23, RPS3, RPS5, 
SFRPI, SM0C2, S0X6, SRPK2, T3-TR-RXR, Ubiquitin, UCP3 

ACATl, Caspase, CD58, CMPKI, ERKI/2, FU5, HMGBI (includes 28 
EG:3I46), lEN Beta, IgG, ILI, ILI2 (complex), IL6ST, ILF3, 
IMPDH2, Insulin, Interferon alpha, Jnk, LDL, LGALSI, LGALS3, 
LRAT, MapK Mek, NFS, P38 MAPK, PCCB, PDGF BB PI3K, PIBFI, 
PRMTI, RABFPl, Ras, SIC25AU, SlAlSa/b, TNIK 

ARHGDIG, ARHGEF ARHGEFl I, CPS I, EMP3, ETEDH, FAU, 27 
FNBPI, GLI3, IFI6, lARPI, LEP Lpa receptor, P2RX7, Pka, 
PD<NB2, PXN, RAGE Ras homolog RCN3, RGNFF RHOH, 
RHOU, RHPNI, RND2, RPL8, RPL26, RPL29 (includes EG:6I59), 
RPLP2, SCAMP2, St3gal, ST3GAL3, TNF UCP3, YWHAZ 

ABHD5, ALDH3A2, BAG3, BCL2, beta-estradlol, BIK, BOG, 25 
C1QTNF6, CDON, CTNNBl, CTSH, CXADR, EEF2, FU, FOLHI, 
GADD45B GIMAP5, H0XC6, HSPB8 IFTI22, IGFBP6, KCNMBI, 
METTL7A, MME MRVII, PLIN2, PUNS, PNPLA2, PTGIS, PTPRK, 
PTPRU, SEMA5A, SFRPI, SULTlAl, SULTIA3 

Arginase, AZGPI, BMP6, C120RFI 1, CI9ORFI0, C2IORF7, 23 
CALM2, C0LI6AI, COL4A6, CSRP2, CTSH, CTSL2, DBNDD2, 
DPYSL3, DYRK2, ENPPI, G0LGA7 HRAS, HTT, IFNG, ILI 3, 
LPHN2, LRBA, LY6E Pdgf, PLODI, RABI2, RAB33A, SLCIA3, 
SRM, SV2A, TCFI2, TGFBI, UNC5B, ZDHHC9 

ADAMIO, AHCYL2, BCL2LM, CABCI, C0R02A, DUT (includes 22 
EG:I854), ERH, Fl IR, GRPR LETMDI, MAPKI, MARCKS (includes 
EG:4082), MESX NFATC2IP NMB NMBR NMTI, NR3CI, PARK7, 
PLK3, PMMI, PRKRA, PRPSAPl, REEP5, RNFI44B, SH3DI9, 
SLC45A2 SMNI, SNUPN, SPOP TP53, TRAF6, UBE2T 
ZNF259, ZNF346 

ABCD3, ACTCI, Actin, ATPI3A5, ATP5J2, ATPase, CDKN2A, 21 
CLPX DDXI9B, DONSON, FAM167A, GPHN, HIPIR HLTF, 
KATNAI, MED27, MIR373, MIR297-Z NCALD, NSF 0RC5L, 
PACRGL, PHACTRi, phosphatidylinosltol 3,4-diphosphate, PLSI, 
P0LR2H, RNA polymerase II, RSPOl, SDCCAG3 (includes 
EG:i0807), SNX16, SUOK TBX3, TERT T0B2, VIM 

ACAA IB, CPT2, CUEDC2, DAP EHHADH, ERBB2, ESRI, FRRSl, 1 9 

GALNTI, GALNT2 GALNT3, GALNT6, HCG 2023776, IGFBP6 
LAMPI, IAMP2, MMSI9, MRCI, MYC NC0A4, 0LFML3, PAICS 
PHF5A, Polypeptide N-acetylgalactosaminyltransferase, PPARG, 
PTRF QKI, RAB34, SCG5, SETDBI, SFRPI, STIMI, TPD52, WWCi 

ALXl, APEH, C20ORFI60, C0PB2, C0X2, C0X4I1, DAGI, EPO, 1 9 

FOS, GMFB GRB2, GRR GSTKI, GYLTLW, HCLSi, HNRNPH2, 
HNRNPR Hydrolase, KHSRR MSTIR ONECUTI, P2RYI, PDE6C 
PDE6G P0P7 RAB3GAPI, RCC2, RNFI28 RPL7, RPS7, RPSW, 
SMARCD2, TMEM62, USP8 YPEL5 



25 



25 



25 



17 



17 



15 



15 



14 



13 



13 
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Table 5 Networks generated from endometrial gene expression data of HF versus LF heifers by IPA (Continued) 



12 Embryonic Development, Tissue Development, ACOTl 3, APBBl, BACEl, C110RF52, CBS, CEtAh Coup-Tf, DCK, 19 13 

Tissue IVIorphology DKKl, DSC2, FOXAl, HMGB2, HNF4A, IER5L, JKAMP, KRRl, LRP, 

LRP2, LRP5, LRP6. MESDC2, NR2F1, NR2F2, PARP4, PHB (in 
dudes EG:5245), PKP2, Plasminogen Activator, P0U5F1, RSPRYl, 
Secretase gamma, 5LCi7A5, TXNDCiZ UBE2D3 
(includes EG:7323), UBE2V1, WNT4 



was indicative of reproductively healthy animals, with 
good heat detection and insemination technique providing 
confidence in retrospective fertility status. However, it is 
important to note, other factors potentially contributing 
to the conception rate differences observed between HF 
and LF heifers, including oocyte quality and oviductal 
environment, were not analysed in this study. 

Endometrial function plays a critical role in pre- 
implantation embryo survival. Consequently, much work 
has focused on the biochemical and molecular phenom- 
ena surrounding the progression of an estrous cycle 
[27]. The present study is novel as it provides informa- 
tion on gene expression during an important period of 
the estrous cycle: the mid-luteal phase, otherwise recog- 
nized as a critical period of embryo loss during pregnancy 
[5,7] between animals of high and low reproductive cap- 
acity. Reiterating the importance of examining transcrip- 
tion during this phase, Salilew-Wondim et al. [31] recently 
found more extensive differential gene expression in endo- 
metrium harvested from heifers on D7 (an estrous cycle 



prior to embryo transfer) between heifers that conceived 
and those that returned to estrus before day 21, when 
compared with D14. 

GALNT6, encoding enzyme UDP-N-acetyl-alpha-D-gal 
actosamine: polypeptide N-acetylgalactosaminyltransferase 
6, was the most abundantly expressed gene in LF heifers. It 
was 6.7 fold up-regulated in LF compared with HF heifers. 
This is the first report of expression of this gene in Bos 
taurus. The GALNT6 gene is located on chromosome 
5 in the bovine genome and shares a coding region with 
SLC4A8, a sodium bicarbonate co-transporter. Expres- 
sion of this gene in humans is implicated in the synthesis 
of oncofetal fibronectin (onfFN) [64], a protein found in 
plasma and cervicovaginal secretions; increased concentra- 
tions of which has been associated with abnormal preg- 
nancy [65,66]. However, Feinberg et al. [67] reported 
increased protein levels of onfFN at the trophoblast-endo- 
metrial ECM interface in human pregnancy tissues from 
gestational day 20 to full term in healthy pregnancies. These 
observations suggest that differential expression of the 
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1.214 




1.385 



1.123 



-1.240 



1.261 



Figure 4 Network #2; cellular growth and proliferation, connective tissue development and function, skeletal and muscular system 
development and function. The network is displayed graphically as nodes (genes). The node color intensity indicates the expression of genes; 
with red representing up-regulation and green, down-regulation in LF versus HF endometrium. The fold value is indicated under each node. 



enzyme GALNT6 may have consequences for embryo sur- 
vival and that this may be time specific however its role is 
currently unclear. 

Pathway analysis is widely used to analyze gene expres- 
sion data and serves as an effective tool for delineating the 
underlying biological processes involved in mRNA aberra- 
tions [68-71]. Biological pathways altered in the current 
study included: cellular growth and proliferation, lipid me- 
tabolism, tissue remodeling, ECM mineralization, inflam- 
mation, angiogenesis, and metabolic exchange. 



Cellular growth and proliferation 

Owing to its regenerative nature, the endometrium under- 
goes highly complex but tightiy regulated cellular prolifera- 
tion and differentiation throughout the estrous cycle [72]. 
There is little published information on the molecular 
mechanism of bovine endometrial proliferation throughout 
the estrous cycle however, studies examining uterine tissue 
of non-pregnant ewes during cycle days 0 to 15 showed an 
increased rate of cellular proliferation between days 0 and 
4, decreasing by day 15, suggesting a proliferative dispos- 
ition is normal earlier in the estrous cycle [73]. Results 
from our study indicate that LF animals could be experien- 
cing an abnormal decline in cellular growth/prolifera- 
tion i.e. 21 genes implicated in cellular proliferation 
inhibition, including FST [74], NPPC [75], GJAl [76], 
S0X6 [77], were up-regulated in the LF animals. Of 
these genes FST, NPPC and GJAl were previously found 



to be expressed in bovine endometrial tissue [78,79]. 
Substantial inhibition of endometrial cellular proliferation 
would retard the development of a secretory endometrium 
and suppress endometrial maturation [80], thus making 
successful implantation unlikely. 



Angiogenesis 

A critical element of tissue growth and development is 
the growth of new blood vessels, also known as angiogen- 
esis [81]. Generally inactive in healthy individuals and ani- 
mals, angiogenesis plays an active role in endometrial 
function, as well as growth of ovarian follicles and CL dur- 
ing the reproductive cycle [82,83]. In a highly proliferating 
tissue such as endometrium, and particularly during 
the hypothesized window of proliferation day 0 to 14/ 
15, angiogenesis is necessary for the provision of nutrients. 
Factors controlling angiogenesis include growth factors, 
nitric oxide and matrix metalloproteinases (MMPs), of 
which MMP19 was down-regulated in the LF animals [84]. 
Also down-regulated, high-mobility group box 1 (HMGBl) 
which codes for a protein which has previously been iden- 
tified in uterine fluid of dairy heifers on day 7 post estrus 
[22] . A role for members of the HMBG family in angiogen- 
esis is supported by their expression during mouse em- 
bryogenesis [85] with higher expression levels found in 
proliferating cells [86] and lower expression in fibroblasts 
from old-age humans [87]. Down-regulation of these and 
other angiogenic genes, which was the case in LF animals, 
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could prevent the necessary angiogenic cascades synergistic 
with cellular proliferation that dominate the mid-luteal 
phase [34]. 

Lipid metabolism 

Lipid metabolism appears in three of the top 5 networks, 
suggesting its importance as a metabolic process in uter- 
ine physiology. Genes involved included ACATl, CCAT, 
LGALSl, PCCB, SRD5A1, FASN and PPARA. In particu- 
lar, increased PPARA transcript abundance, as observed 
in LF heifers, coincides with increased fatty acid catabol- 
ism [88]. Fatty acids are essential precursors to steroids 
and eicosanoids, metabolites necessary for normal ovar- 
ian and uterine function [9]. Furthermore, studies have 
shown fatty acid supplementation positively influences 
reproductive performance [9,54]. 

Fatty acid synthase (FASN) exhibits its anabolic cap- 
acity by aiding in the conversion of dietary carbohydrate 
to fat, which is subsequently organized into hepatic adi- 
pocytes and lactating mammary tissue as triglyceride and 
milk lipids, respectively [89,90]. It has also been found that 
expression of FASN peaks during the proliferative phases 
of the menstrual cycle [91]. Metabolic demands are par- 
ticularly high during this phase as a result of the exten- 
sive endometrial remodeling and reconstruction, a central 
theme to both the estrous and menstrual cycles. Increased 
FASN would be favorable in such a demanding situation to 
deliver the required fatty acid for the assembly of new cell 
membranes, modification of DNA transcriptional machin- 
ery and hormone construction. Interestingly, expression of 
FASN was down-regulated in the LF heifers suggesting the 
aforementioned processes were compromised in these ani- 
mals, which potentially affecting their ability to conceive. 

Steroid 5a-reductase type 1 enzyme is involved in the 
metabolism of progesterone that is found in uterine and 
cervical cavities. Murine gene knock-out studies have 
shown that parturition is adversely affected by aberrant 
expression of this gene, impeding cervical ripening and 
fetal delivery as a result of elevated progesterone levels 
in the cervix [92]. Expression of the gene coding for this 
enzyme was up-regulated in LF heifers, thus progester- 
one catabolism is likely to be active in these animals. As 
high progesterone levels are positively associated with 
embryo survival [60,93], it is therefore possible that the 
LF animals are experiencing low local progesterone con- 
centrations and ultimately, this could be contributing to 
their low conception rates. 

Cellular and tissue morphology and development 

The ability of cells to generate alternate cell types whose 
phenotype is different from that of the source tissue is 
known as plasticity. Endometrial epithelial and stromal cell 
proliferation, as discussed previously, is a complex multi- 
component process involving cues from extra-cellular 



growth factors and ovarian hormones [72,94]. However, in 
their absence, isolated bovine endometrial stromal cells 
exhibit the ability to develop into bone [95]. Results from 
our microarray study showed a large representation from 
this biological category, with 36 DEG enriched. Genes im- 
plicated in cell and tissue morphology and development 
which were down-regulated in low fertility heifers in- 
cluded, PPARA, IL6ST, GJAl, SFRPl and IL-33. 

One particular biochemical pathway which facilitates 
cellular transformation includes extracellular matrix 
mineralization (ECM) [96]. A well known regulator of 
ECM mineralization is the activin a-FST system. Activin 
A inhibits ECM mineralization whereas FST, an activin 
antagonist which prevents activin-receptor interaction 
[97], increases mineralization in cell cultures [98]. Trans- 
genic female mice with gain-of-function FST, in which 
mouse foUistatin was over-expressed, developed thin uteri 
and small ovaries, resulting in infertility [99]. FST was dif- 
ferentially expressed between HF and LF heifers, indicat- 
ing a role for this gene pathway in mid-luteal endometrial 
homeostasis and early embryo survival. 

ECM remodeling, occurring during both pregnancy and 
the estrous cycle, facilitated by the matrix-metalloprotei 
nases, ensures the provision of a suitable structural micro- 
environment where the embryo can grow [100,101]. 
Matrix-metalloproteinase-19 {MMP-19), an important 
molecule in this pathway and which was down-regulated 
in LF heifers, plays a significant role in ECM remodelling 
[102]. Interestingly, Wathes et al. [103] reported that dif- 
ferential expression of genes MMP - 1, 2, 3, and 13 two 
week post partum in the bovine endometrium, was highly 
correlated with differential expression of IGF binding pro- 
tein 4, a known antagonist of IGFl expression [104]. The 
IGF system, in particular IGFl, is associated with several 
reproductive processes in cattle including preimplantation 
embryo development [105-107]. 

The transforming growth factor |3s (TGF-|3) are multi- 
functional cytokines that also regulate tissue remodelling 
and repair [108,109]. High expression of TGF-|3 has been 
observed during pro-estrus and diestrus [110,111] thereby 
highlighting the role for TGF-ps in endometrial remodel- 
ling, an important process impeding estrous cycle transi- 
tion [112]. Transforming growth factor beta 1 induced 
transcript (TGFBIII) was down-regulated in the LF ani- 
mals, suggesting altered or irregular endometrial remodel- 
ling in these animals which may be contributing to the 
conception rate differences observed between the two di- 
vergent fertility groups. 

Inflammation 

Inflammation is an innate cyclical physiological process 
facilitating progression of reproductive cycles in the endo- 
metrium. The animal model in this study isparticularly 
useful for the identification of inflammatory pathways 
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associated with uterine low fertility for numerous reasons. 
Firstly, there has been no mitogenic challenge. This study 
strictly examines gene expression between high and low 
conception rate animals without influence from any ex- 
ogenous metabolites, either dietary or pharmaceutical. Sec- 
ondly, tissue sampling occurred during an estrous cycle 
where no embryo was present. Lastly, the study employed 
nuUiparous heifers where the likelihood of uterine infection 
is low, as was demonstrated by the lack of clinical evidence 
of metritis, endometritis, pyometra or metaplasia across all 
heifers examined. 

In total 55 DEG featured in inflammatory linked path- 
ways. It is clear from the high proportion of DEG that in- 
flammation is a central theme in estrous cycle and uterine 
sub-fertility physiology. IL-33, a cytokine which influences 
the production of other pro-inflammatory cytokines IL-5, 
IL-13 and chemokine GM-CSF [113] was more highly 
expressed in LP animals. In addition, IL-33 regulates tran- 
scription of endothelial cells in inflamed rheumatic tissues 
[114]. As mentioned previously, cell plasticity is altered in 
a state of chronic inflammation or trauma. Inflammation 
due to up-regulated IL-33 could be altering the constitu- 
tion of the endometrium in the LF animals, and thus im- 
peding embryo implantation. Hence, low conception rates 
could be directly linked to inflammation induced, altered 
cellular plasticity, in uterine endometrial tissues. 

Metabolic exchange 

Similar to Forde et al. [33], Bauersachs et al. [29] and 
Salilew-Wondim et al. [31], genes coding metabolite trans- 
porters, specifically the solute carrier (SLC) family mem- 
bers were found to be differentially expressed between HF 
and LF animals. The five SLC genes identified were; 
SLC1A3, SLC17AS, SLC25A12, SLC25A24, SLC4SA2. The 
most abundantly expressed gene of the entire DEG list, 
SLC4SA2, was 8-fold more highly expressed in the uterus 
of LF relative to HF heifers. As the name suggests SLC 
genes are involved in the transfer of solutes across the 
cell membrane, particularly amino acids [115-117]. Amino 
acids are fundamental for the normal growth and develop- 
ment of the early embryo, acting as precursors of nucleic 
acids and proteins, osmolytes and signaling molecules. 
Concentrations of amino acids in oviductal and uterine 
fluid during the estrous cycle have been reported to modu- 
late with stage of cycle, systemic progesterone environment 
and differ compared with plasma, demonstrating their ac- 
tive transport in these tissues [118-121]. The endometrium 
functions as a secretory layer, suggesting the importance of 
metabolite exchange in this specific tissue. Animals with 
less efficient metabolic exchange in the uterus may be 
unable to sustain embryo development during early 
pregnancy, and thus be experiencing recurring early 
embryo loss. 



Microarray analysis was carried out on endometrial 
tissue, an amalgam of varying cell types. Examining tis- 
sue mRNA gene expression provides an insight into the 
genetic regulation of multiple cell types from the host. It 
was essential to use RNA from all endometrial cell types 
as it is not apparent, as of yet, whether or which individ- 
ual endometrial cell types are contributing to low con- 
ception rates in cattle. Investigations into the types and 
locations of contributing cell types via in situ hybridisa- 
tion or immunofluorescence would assist in the develop- 
ment of proposed hypotheses. 

Conclusion 

Global endometrial gene expression profiles during the 
mid-luteal phase of the estrous cycle, in HF and LF heifers 
was investigated, and the most significant biological path- 
ways likely to be involved in uterine function and embryo 
survival identified. The new knowledge generated offers 
substantial insight into some of the molecular mecha- 
nisms underlying uterine endometrial function and uter- 
ine mediated low-fertility, during the early to mid luteal 
phase of the estrous cycle in cattle. Furthermore, expres- 
sion analysis provides invaluable data on key differentially 
expressed genes which may be selected for future SNP 
discovery analysis which following validation may be used 
as genetic markers for fertility and incorporated into 
breeding programmes. 
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